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Abstract 

Supersymmetric Gauge-Higgs Unification is a well-motivated new physics scenario, 
both in heterotic model building and from the perspective of higher-dimensional 
Grand Unified Theories. When combined with radion mediated supersymmetry 
breaking, it allows for very specific predictions concerning the high-scale parameters 
of the MSSM. Using an appropriately modified version of a standard RGE evolution 
code (SuSpect), we derive low-scale predictions which can be tested at the LHC. The 
phenomenological success of our setting depends crucially on the 5d Chern-Simons 
term, which has not been used in previous, less encouraging studies of gauge-Higgs 
unification in supersymmetry. 
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1 Introduction 



If the LHC discovers supersymmetry, it will be a great challenge to relate the measured 
electroweak-scale parameters of the supersymmetric Standard Model to more fundamental 
structures at a high energy scale. One of the best-motivated physics proposals in this 
context is supersymmetric grand unification [1,2]. However, SUSY grand unification by 
itself places only rather limited constraints on the low-energy parameter space. 

An elegant and natural further idea making this setting more predictive is gauge- 
Higgs unification (GHU). In models of this type, some or all of the Higgs scalars are 
extra-dimensional components of gauge fields GHU is for instance realized in many grand- 
unified models derived from heterotic string theory, where one or both of the MSSM Higgs 
doublets can come from the untwisted sector (see e.g. [4]; for a recent review see [5]). 
At a simpler level, GHU can be realized in purely field-theoretic 5d or 6d orbifold GUT 
models [6]. These can be viewed as effective unified field theories, valid directly below 
the heterotic string scale. Such constructions receive independent support from the string- 
scale/GUT-scale problem [7] as follows: One of the possibilities for overcoming this problem 
is the compactification on anisotropic orbifolds [8-11], where one or two of the compact- 
ification radii are much larger than the string length scale. This naturally allows for an 
intermediate effective description in terms of a 5d or 6d orbifold GUT. 

The present paper, which is mainly phenomenologically oriented, does not depend on 
specific string-theoretic realizations of GHU. Our analysis relies only on simple SUSY field 
theory models in which GHU can arise, and whose 4d low-energy effective field theory 
is the MSSM. The earliest and simplest construction of this type is the 5d SU(6) model 
of Burdman and Nomura [12], which we will largely follow. We expect, however, that our 
phenomenological results will carry over to similar models, including more elaborate string- 
derived constructions. Related models include, e.g., the 5d SU(6) model with warped extra 
dimensions of [13] and the 6d models of [14-16]. For more work on SUSY GHU in orbifold 
GUTs, see for instance [17] and references therein. 

The main point of 5d GHU models relevant to our work is easily explained: 5d gauge 
symmetry enforces a Kahler potential in the Higgs sector of the form [18] (see also [19-21]) 

j d"^ j ^^yH{T,T) (Hi + H2){H2 + Hi) . (1) 

Here Hi and H2 are the MSSM Higgs superfields, arising from the adjoint of the 5d gauge 
group after its breaking to the Standard Model group. The prefactor Yh is a real function 
of the radion superfield T, and ip is the chiral compensator of 4d supergravity. If T and 
Lp develop F-term VEVs [22], Eq. ([T]) clearly induces a supersymmetric jj, term as well as 
soft Higgs mass terms. They satisfy the relation 

2 2 I 2| rn\ 

■nil — ^2 — l^^sl 

at the GUT scale, where 2 = + "^l^i 2 diagonal entries of the Higgs mass 

matrix and m\ = B/j, is the off-diagonal element. This is a specific realization of the 
Giudice-Masiero mechanism [23] or of its string-theoretic version [24]. 

^In fact, GHU has a long history in non-supersymmetric models without grand unification (see e.g. [3] 
and many subsequent papers). Here, we take the SUSY GUT idea as our main paradigm. 
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At the TeV scale, the famihar conditions for electroweak symmetry breaking and vac- 
uum stabihty read: 

m\ml - < 0, 

Renormahzation group (RG) running must thus turn the equahties of ^ at the UV scale 
into the inequalities of ([s]) at the weak scale. A numerical analysis is required in order to 
find out whether this is possible at all and, if so, whether it is possible within a realistic 
model. This analysis has to take into account some additional predictions of 5d GHU 
models. In particular, there are strict relations between the Higgs mass parameters, the 
gaugino mass, and the dominant soft sfermion masses and trilinear couplings. 

If the function comes entirely from the gauge-kinetic term of the 5d super Yang- 
Mills action, one finds Yh ~ 1/(T + T). In this case, the relations between the soft pa- 
rameters turn out to be rather restrictive and a realistic low-energy spectrum cannot be 
obtained without extreme fine-tuning [18]. 

If we also include a supersymmetric 5d Chern-Simons term and allow for a VEV of 
the chiral adjoint in the 5d gauge multiplet, an extra contribution ~ 1/(T + T)^ to the 
function Yh{T,T) arises [25]. The Chern-Simons term is generically present in a 5d super 
Yang-Mills theory compactified on S^/7j2, and its coefficient is determined by anomaly 
cancellation [26,27]. It depends on the full field content of the model, which we do not 
fix completely. In particular, the fields of the two light generations can be distributed in 
various ways between bulk and branes. The prefactor of the 1/(T + T) ^-contribution to 
Yh depends on the coefficient of the Chern-Simons term and on the size of the adjoint 
VEV, which is also unknown. Hence we treat this prefactor as an extra parameter. 

The main point of the present paper is to demonstrate that, allowing for this Chern- 
Simons term contribution to Eq. ([T|, a completely realistic low-energy phenomenology can 
be obtained. This scenario is rather constrained since all MSSM soft parameters are given 
in terms of the VEVs of F'^ and F''', a dimensionless parameter c' characterizing the effect 
of the Chern-Simons term, and two mixing angles related to the 5d-origin of the third- 
generation quarks and leptons. We analyze the low-energy phenomenology of this setting 
and discuss observational consequences for the LHC. 

To this end, we numerically solve the renormahzation group evolution of the MSSM 
parameters, with GUT-scale boundary conditions as determined by our GHU model. 
This procedure has been implemented within the public MSSM spectrum generator code 
Suspect 2.41 [28]. We perform two extensive parameter scans. The first one uses a rather 
crude estimate for boundary conditions in the sfermion sector. This is essentially a gen- 
eralization of the analysis of [18], now including the effects of a Chern-Simons term and 
using a state-of-the-art RG code. The second scan uses realistic boundary conditions, de- 
rived from the Burdman-Nomura model [12]. It is somewhat more involved because of the 
specific relations between 4d Yukawa couplings and fundamental model parameters. 

In both cases we find regions of parameter space where all present experimental bounds 
are satisfied. In particular, we can have sufficiently high Higgs and sparticle masses to evade 
direct search bounds, comply with i?-physics constraints on rare decays, and also obtain a 
neutralino as the lightest supersymmetric particle (LSP) with a dark matter relic density 
compatible with WMAP results. We conclude that SUSY gauge-Higgs unification with 
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radion mediated SUSY breaking is indeed a phenomenologically viable scenario. We also 
point out some characteristic LHC signatures which seem to be typical for this framework. 

The paper is organized as follows: In Section [2| we explain the appearance of the GHU 
relation, Eq. (|2]), in a large class of models. In Section [3| we present a concrete 5d model 
and give the formulas for the GUT-scale gauge- and Higgs-sector soft terms as functions 
of the model parameters. Section |4] contains a discussion of the expected RG running 
behaviour, particularly of m^, and of its consequences for identifying phenomenologically 
promising parameter space regions. We discuss our general setup for numerically analyzing 
the RGEs in Section [5j Assuming a simplified set of sfermion soft terms, we present a 
first such analysis in Section [6j We proceed in Section [7] by explaining how a fully realistic 
sfermion sector can be included in the analysis. The corresponding numerical results are 
presented in Section |8] and Conclusions are given in Section [9j 



2 The origin of the GHU relation for Higgs mass pa- 
rameters 

In string-derived or orbifold GUT models, the MSSM Lagrangian will generically depend 
on several moduli fields. We focus on models where SUSY breaking is communicated to 
the MSSM predominantly through the moduli. Furthermore, we assume that the Kahler 
potential depends on the Higgs superfields Hi and H2 only in the combination Hi + H2 
and Hi + H2- This assumption will be justified momentarily for a certain class of models. 
The leading part of the Higgs action then reads 

Sd [ d^x [ d^e Yh{Z', Z^) (Hi + H2)(H2 + Hi) , (4) 



where Yh is some real analytic function of the moduli fields . 

In GHU models there is no superpotential contribution to the n term. The Higgs mass 
parameters are therefore entirely determined by Eq. Q: For canonically normalized Higgs 
fields they are given by 

rriH, = = -F'F'didjlogYH, 

±fx = F^ + F'djlogYH, (5) 
±Bi2 = + F^di hgYuf - F^F'^didjlogYn. 

These equations obviously imply Eq. 

If we start from a 5d model, it is straightforward to see why the Kahler potential 
always depends only on the combination Hi + H2 (and its complex conjugate): Recall that 



^ Note the sign ambiguity in n and -B/x: Simultaneously changing the signs of both fi and corresponds 
to a redefinition of one of the Higgs fields, say Hi, by Hi —Hi. The signs of the Yukawa couplings 
can be kept unchanged by an analogous redefinition of the right-handed matter fields it couples to. The 
overall sign of ^ and B/i cancels in the RGEs. Therefore, given a certain high-scale model, a simultaneous 
sign change in the last two lines of Eq. ^ will leave the absolute values of all weak-scale parameters 
unchanged, flipping only the signs of fi and Bfi at the weak scale. The convention that Bfi is positive at 
the weak scale then determines the signs in Eq. (jsj). 
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a generic 5d super- Yang-Mills theory contains, in terms of 4d superfields [29-31], a vector 
superfield and a chiral adjoint $. The fifth component of the gauge field, A^, forms the 
imaginary part of the scalar component of In the 5d kinetic action, the chiral adjoint 
appears only in the combination $ + $^ [31]. The reason is that, in this combination, 
drops out of the lowest component of the real superfield $ + "l*^, ensuring the absence of 
non-derivative couplings of A^^ 

The decomposition of the adjoint of the 5d gauge group G in irreducible Standard 
Model representations is vector-like and can contain pairs of weak doublets: 

Ad(G) ^(l,2)_i/2 + (l, 2)1/2 + ... (6) 

If this is the case, and if the zero-modes of these doublets in $ survive compactification 
and GUT symmetry breaking, they can be identified with the Higgs fields Hi and H2. It 
is then only the combination Hi + H2 which appears in the 4d Kahler potential. 

Intriguingly, this peculiar combination of chiral superfields has also been found in 
heterotic Eg x Eg orbifold models, where no use of an intermediate 5d effective theory is 
made [19-21]. We briefiy describe the situation following [20]: The Kahler potential can be 
expanded as 

K = IC + yaaA^A"" + yp-pB^B~^ + [Z^pA'^B^ + h.c.) + . . . (7) 

Here /C, 3^, Z are functions of the moduli, and A" and B^ are matter fields transforming 
in the 27 and 27 of Eq, respectively. Suppose that the MSSM Higgs fields Hi and H2 are 
contained in two such fields A and B. For the desired combination Hi + H2 to appear, one 
needs 

yAA = ysB = Zab = Zxb ■ (8) 

Indeed, for untwisted matter fields A, B associated with a common complex plane, it was 
found that ^ 

yAA = yBB = Zab = Zab = ^t){U + U) ' 

where U and T are the complex structure and Kahler modulus of that plane. By contrast, 
for twisted matter fields, or for untwisted matter fields associated with distinct planes, or 
with a plane without U modulus, one has Zab = 0. 

We now argue that the above string-theoretic results are closely related to the previous 
5d argument. Firstly, if A and B are associated with the same complex plane, one can take 
a 6d limit in which the two corresponding compactification radii remain large. In this limit, 
the fields A and B are the extra-dimensional components of the 6d gauge field. Secondly, 
the presence of a shape modulus U of the corresponding large allow for 5d limit, in 
which the compact space becomes an interval. Thus, the string-theoretic conditions which 
ensure that the combination Hi + H2 appears are precisely those which are needed for an 
appropriate 5d limit to exist. Even in regions of moduli space which do not correspond 
to that limit, the structure enforced by 5d gauge invariance survives. In other words, the 
simple 5d argument given earlier appears to be sufficient to understand the situation in 



■^In fact, an analogous argument forces the Kahler potential to depend on the radion modulus T only 
through the combination T + T. In this case, the imaginary part is the 5th component of the graviphoton 
of 5d supergravity, and the real combination T + T ensures the absence of non-derivative couplings of this 
gauge-field component [29]. 
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heterotic orbifold models as well. It would be interesting to work this out in more detail, 
which is however beyond the scope of the present, mainly phenomenologically oriented 
paper. 

As we have already emphasized in the Introduction, and as will become clear in sub- 
sequent sections, the 5d supersymmetric Chern-Simons term plays a central role in our 
analysis. Such a term cannot arise in the tree-level dimensional reduction of a theory with 
more than five dimensions. However, it is generically induced at one loop [32]. More specif- 
ically, the radiative generation of a Chern-Simons term in compactifications from 6d to 5d 
was discussed in [33]. Thus, the 5d Chern-Simons term is consistent with a lOd origin of 
the theory. 



3 A concrete 5d realization 

We now turn to the phenomenological prospects of concrete 5d models. Our main example 
will be a generalization of the SU(6) model of Burdman and Nomura [12], with a larger 
gauge group containing at least U(6) = SU(6) xU(l). The 5d gauge theory is compactified 
on /{'L2 X Z2). The only relevant modulus is the radion superfield T = p + iB^ + . . . , 
where is the fifth component of the graviphoton. The real part p is normalized such 
that 2'Kp is the volume of the original S^. We assume that it is eventually stabilized at 
(p) = R. The Higgs doublets are contained in the superfield $ = S + iA^ + . . . (where S 
is the chiral adjoint of the 5d gauge multiplet, and is the fifth component of the gauge 
boson). 

The U(l) gauge factor is assumed to be broken on the boundary. Furthermore, orbifold 
boundary conditions for the gauge fields can be chosen such [12] that the remaining SU(6) 
is broken in 4d to the Standard Modelj^ and that the only components of $ with zero 
modes are the Higgs fields. Their origin is particularly obvious in the SU(5)-decomposition 

35 = 24 + 5 + 5 + 1 (10) 

of the adjoint of SU(6). Our fields Hi and H2 are the doublets contained in the 5 and 5 
respectively. 

Before orbifolding, the 4d effective theory contains the full gauge multiplet and chiral 
adjoint $ as well as the radion T. The corresponding leading-order action, which has been 
analyzed in [29], contains a term 

^^S/^'" /^'^^^^tr($ + $t)^ (11) 

Retaining only the $ components that survive the orbifold projection, i.e. the Higgs fields, 
this becomes [18] 

SD^-^Jd'xJ d'eW7^(Hi + H,){H, + Hi). (12) 

Note that the coupling of T to the Higgs fields depends on the choice of Kahler-Weyl 
frame. We work in a frame where K = — 31og(T + T + . . .). The 4d metric has not been 
rescaled after compactification. 

Apart from an extra U(l) which is Higgsed by a brane field. 
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The 5d theory will in general also contain a Chern-Simons term. Its supersymmetrized 
version includes a cubic term in $, which couples to the radion according to [25] 

S^'-^jd-xj d^e^^ (^)^^($ + $)^ (13) 

After orbifolding and allowing for a non-zero expectation value ($) = v Ig, this part of the 
Chern-Simons action contributes as 

SD^-^Jd'xJ d'ew {-^J + + (14) 

to the quadratic Higgs Lagrangianj^ Here we have introduced the dimensionless constant 
c' = 2cvg'^. Note that the group U(6) = SU(6) xU(l), which we use here and below, is only 
the simplest extension of SU(6) allowing for an SU(6)-preserving $-VEV|^ Larger and, in 
particular, simple groups, such as SU(7), are clearly possible. 

We regard C clS db free parameter for the purpose of our analysis: While boundary- 
anomaly cancellation fixes c for a given field content, we specify neither this field content 
nor the value of v. In particular, different distributions of the light generations between 
the two branes affect the value of c. However, positivity of the kinetic terms (given just 
below) leads to the constraint 

c>-l. (15) 

The MSSM scalar potential for canonically normalized Higgs fields (which we also call 
Hi and H2 by abuse of notation) reads, at quadratic order, 

V = m?|ifi|2 + m^li/aT + ml{H2Hi + h.c.) . (16) 

The mass parameters can now be calculated from the Higgs kinetic function 

y„(T.T)^^(l + A)^. (17) 



gi \ T + TJ T + T 



which follows from Eqs. (12) and (14). According to Eqs. ([s]), the parameters [i and 
read [25] 

f'^ 1 + 2c- 

^-^ = ^ -2i?TT^' 

2 2 2 ,™,2 (F^F^ + h.c.) 1 + 2c' IF^P 2c'^ , , 

^ 2 3 I I 1 + c' (2i?)2 (1 + d f ' ^ ^ 

where we have introduced a parameter e^^ = ±1 to account for the sign ambiguity in 
m| = B\i and [i. 

From the 5d gauge-kinetic and Chern-Simons action we obtain, after dimensional 
reduction, 

S j d^x j d'e'^ tr WW^ + h.c. + 2c7rR J d^x j d'O tr ($ WW^) + h.c. , (20) 

^Our <I>-VEV V — V4d is a VEV in the 4d effective theory, ft is related to the corresponding VEV 
of the underlying 5d theory by v^d — {p/R)'>^5d- Here the mass dimensions of the 4d and 5d scalar fields 
are the same since we assume that the complete leading-order 5d Lagrangian has a prefactor . 

^The significantly more complicated case of SU(6)-breaking expectation values is considered in [25]. 
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which eventually gives the 4d gauge-kinetic term 

Sd'^ j d^x j d^e + c'^ tiWWo, + h.c. . (21) 
It determines the gaugino masses 



-T 

F 1 



as well as the 4d gauge coupling 



Ml/2 = (22) 



.2 + (23) 
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The soft masses and trilinear terms for the matter multiplets are more model- 
dependent. Quite generally the relevant piece of the kinetic action can be written as 

[d'x [ d'9^^\Yu{T,T)\U\^ + Yq{T,T)\Q\^ + Yn{T,T)\D\' 
J J ^ (24) 

+ Ye{T,T) |Ep + Fi(T, T)|L|2 + Yn{T, T)\N\'' . 

The kinetic functions Yx (with X = U,D, Q, E, N, L standing for up-type and down-type 
right-handed quarks, quark doublets, charged and uncharged right-handed leptons and 
lepton doublets) determine the soft masses according to 

r)2 

ml = -\F^\^—=\ogYxiT,T). (25) 
The trilinear couplings are given by 

Au,D = F^-^ log {YhYqYu,d) , (26) 

AE = F^-^\ogiYHYLYE). (27) 
Note that we define the A term with a negative sign in the Lagrangian: 

CD- [AuyuH2QU + ADyoHiQi) + AeVeHiLE + h.c.) . (28) 

The precise form of the matter kinetic functions depends on the model under con- 
sideration. We will assume throughout that the first two generations of MSSM matter 
are brane-localized and that their GUT-scale soft terms are negligible. This gives no-scale 
boundary conditions for the first and second generation. For the third generation, we will 
consider two cases: first the approximation that only the top quark receives a Yukawa 
coupling induced by the 5d gauge coupling, and second the case of the Burdman-Nomura 
model with realistic top, bottom and tau Yukawa couplings. 



^ Note that this agrees with [25] only after a substitution c' c'/2, which is due to our modified 



definition of c'. However, after this substitution, it becomes apparent that Eqs. (18 1 and (19 1 are truly 
different from [25]. This is the result of our SU(6)-preserving <i>-VEV, as opposed to the SU(6)-breaking 
$-VEV of [25]. 
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4 Expected running patterns 



Before we present the numerical results, let us briefly discuss what general features we 
expect. Needless to say, the complete system of two-loop renormalization group equations 
(RGEs), which will be solved numerically in the following sections, is far too complicated to 
permit an analytical treatment. Some aspects can nevertheless be qualitatively understood 
by inspection of the dominant contributions to the one-loop RGEs. 

The scale of electroweak symmetry breaking (EWSB) is defined as usual as the geo- 
metric mean of the stop masses, Mewsb = ^mj^m^^- At this scale the conditions of Eq. ^ 
have to hold. Once EWSB occurs, one finds the well-known relations between /i, Bfi = m|, 
Mz, the Higgs soft masses mj^, and the ratio tan/? of Higgs expectation values: 

/i^ = - [tan 2/? (m^2 tan /? - m^^ cot p) - M|] , 

I (29) 
5/i = - sin 2/3 [mjj^ + mjj^ + 2/^^] . 



We focus on the region of moderately large tan (3, roughly tan /3 > 5, to ensure that 
the tree- level bound on the lightest Higgs mass, nih < Mz, is approximately saturated. 
The Higgs mass can then be lifted above the direct search limit by radiative corrections, 
mainly due to stop loops. 

The latter involves a significant fine-tuning (the notorious MSSM "little hierarchy 
problem" ) , because the soft mass scale must be large compared to Mz instead of being of 
the same order of magnitude, which would be the natural situation. For sizeable tan (3 one 
has 

Mi ^ 

-Y ~ -^2, (30) 

so must be negative and small compared to typical soft masses. We will not discuss 
this fine-tuning any further (see however [34]), but accept it and focus on the implications 
for models with GHU boundary conditions. One immediate consequence is that vn^ > 
at the GUT scale, because the soft mass and hence also runs down towards lower 
energies (the running of fi is insignificant). While this also fixes the GUT-scale sign of 
mf to be positive, either sign for Bjj, is possible (cf. Eq. |5|. In other words, we can have 
ej^ = -|-1 or e// = — 1 in the GHU relations 

<=< = %i?/^-H^ (31) 



However, as we will now argue, en is always determined by the sign of jj: Out of the 
four sign choices > or /i < and e^j^ = ±1, only two can generically lead to realistic 
spectra. To establish this observe first that ml will typically not evolve by more than a 
factor of 2 — 3, and therefore remains of the order of magnitude of the typical soft mass 
scale during KG running. Furthermore, we just stated that m"^ at Mewsb should be small 
compared to the typical soft mass scale, and that tan/5 should at least be moderately 
large. From all this it follows that S/x at the EWSB scale should be small compared to the 
typical soft mass-squared scale as well, as can be read off from 

ml + , , 

tan p + cot(3= ^— — ^ (32) 
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(which is equivalent to the second hne of Eqs. (29 )). We will now show that requiring small 



EWSB-scale Bfi generically fixes ejj in terms of sign(/i). 

The RG evolution of is primarily governed by the terms involving the top trilinear 
couphng and the weak gaugino mass: 

167r'^S/i = fi{6At\yt\^ + QglM^) + ... (33) 
We can choose positive gaugino masses without loss of generality. Let us now discuss the 



relevance of the two dominant terms on the r.h. side of Eq. (33): 

The gluino contribution to the At-KGE forces At to run negative towards low scales. 
This is fairly universal, i.e. more or less independent of the values of the other parameters. 
The value of At at any given scale is thus to a good approximation dictated only by its 
GUT-scale boundary value and Mi/2- The running of gaugino masses is also approximately 
universal: at one-loop, they simply evolve according to the respective gauge coupling beta 
functions. In the RG evolution of S/z, At will therefore always dominate at low energies, 
when it has become large and negative and when also yt has grown large. Correspondingly, 



the M2 term on the r.h. side of Eq. (33) can dominate only at energies near Mgut, before 
it is overwhelmed by At. 

For negative yU, initially increases from its GUT-scale value and then runs down; 
for positive /z, it evolves in the opposite way. The relative importance of the At and the M2 
contributions is set by their GUT-scale initial values: the larger At at Mqut, the longer 
it will take to run negative and to finally dominate the RG evolution. For small or 
negative GUT-scale At, the B^x running at low energies is more important than the initial, 
M2-dominated phase near Mgut- 

The direction and slope of the running of B^ are set by the sign and magnitude of /x, 
which itself does not run significantly as mentioned. We observed before that Bn should be 
small at the EWSB scale — for the sake of the argument, let us try to construct a situation 
where it is exactly zero. It should in particular change significantly with respect to its initial 
GUT-scale value, so should be sizeable. Furthermore, changing the sign of /i will lead to 
B^ evolving in the opposite way (at least as far as the evolution is governed by the terms in 



Eq. (33)). If there is a solution with, e.g., en = —1 and sign(/i) = —1, we therefore expect 
a nearby mirror solution for the opposite sign choice. On the other hand, changing only 
one of the signs will generically not lead to a solution due to the approximately universal 
behaviour of At and M2. We have sketched this behaviour in Fig. [T]for large GUT-scale 
At, and in Fig. [2] for small or negative GUT-scale At. 

Thus, while one might naively have expected four branches of solutions of the RGEs 
to give realistic spectra (corresponding to the two choices of each sign(/i) and of e^^), by 
the above discussion there should in fact appear only two. Furthermore, we expect sizeable 

in all cases. 
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Figure 1: Qualitative RG evolution of ml (blue dot-dashed curve), (red dashed curve), 
and ml (green solid curve) as a function of the scale s, between s = Mewsb and s = Mqut- 
At at the GUT scale is sizeable and positive. The green dotted curve is m^ for the wrong 
sign(/i), which does not lead to realistic EWSB. Left panel: en = —1 requires sign(/i) = — 1. 
Right panel: en = +1 requires sign(/i) = +1. 




Figure 2: Same as Fig. [T]for small or negative GUT-scale At. Left panel: en = —1 requires 
sign(/i) = +1. Right panel: en = +1 requires sign(/x) = —1. 



5 Numerical analysis: General setup 



For the numerical analysis, we make use of the public state-of-the-art SUSY spectrum code 
Suspect 2.41 [28], appropriately modified to be applied to our SUSY GHU model. The 
usual procedure in Suspect and in other SUSY spectrum codes is to use tan/3 and Mz 



as inputs and to compute fi and 5/i from the EWSB condition Eqs. (29). In our model 



however, the Higgs soft masses and /i and B/j are not independent, since they are related 



by the GHU conditions (31) at the GUT scale. 



The free parameters in the gauge-Higgs sector of the model are actually F'^ /2R, F'^ and 
c', from which the GUT-scale values for M1/2, /i, Bfi and rn'jj_^^ are determined according 
to Eqs. (18), (19) and (22). Together with the GUT-scale values for the sfermion mass 



parameters and trilinear couplings, they furnish a set of GUT-scale boundary conditions 
for the MSSM renormalization group equations. 



It is in principle possible to change the usual procedure of spectrum computation such 
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that /X and become high-scale inputs, while tan (3 as well as Mz are output determined 



by Eq. (29). We have implemented this scheme in Suspect 2.41; the requirement to find 
the correct experimental value of Mz, however, makes parameter scans very inefficient. 

For the present analysis we have therefore chosen a different approach: We work with 
the conventional SUGRA scheme of Suspect 2.41, which takes ta.n l3{Mz) together with 
the GUT-scale values of all soft-breaking parameters except Bfi as input. We only modify 
this scheme by not specifying fixed GUT-scale values for the Higgs soft masses m'jj_^ and 



mfj^, but instead determining them from the GHU boundary conditions Eq. (31). 

Our input parameters are thus Mi/2(Mgut) and tan/5(M^), the two sign coefficients 
sign(/i) and e^j, plus the sfermion mass parameters and A-terms at Mqut- The values of /i. 



Bfi, mjj_^, rn\j^ are computed iteratively applying Eqs. (29) at the EWSB scale and Eq. ( [3l| 
and the GUT scale. When a stable solution is found, the model parameters F^/2i?, F'f 
and d are inferred from M1/2, A* and at Mqut by inverting Eqs. (18), (19) and (22). 



A complication arises, however, from the sfermion sector. As discussed in Section [3} we 
assume no-scale boundary conditions, i.e. a common scalar mass mo = and a common 
trilinear coupling Aq = 0, for squarks and sleptons of the first two generations. The soft 
terms of the third generation, on the other hand, can be non-zero. To be more precise, 
they will depend on F'^ /2R and c' (and possibly also on other model parameters) accord- 
ing to their kinetic functions. This requires an extra level of iteration in the spectrum 
computation. 

It turns out to be convenient to let this iteration act on c'. We thus start the procedure 
described above with an initial guess of c', which is kept constant until a first convergence 
of the spectrum is reached. This has the virtue that the GUT-scale sfermion soft masses 
are unambiguously fixed in terms of c', M1/2, and other input parameters (as will become 
clear once we describe how we are modelling the matter sector) so the EWSB scale does 
not change too much in each iteration step, which could lead to numerical instabilities. 
When convergence is reached, an updated value of c' as computed from M1/2, /i and S/i 
is taken as the new input c', and the whole procedure is iterated until c' converges as well. 

Let us finally list the Standard Model (SM) input values and experimental constraints. 
For the SM input values, we take a-^(Mz) = 127.934, a,{Mz) = 0.1172 and mfc(mfe) = 
4.25 GeV in the MS scheme, and an onshell top mass of rrit = 172.4 GeV [35]. Moreover, 
Mz = 91.187 and = 1.777 GeV, and Gp = 1.16639 ■ 10"^ GeV'^. 

To take into account the limits from direct SUSY searches at LEP [36], we require 
m-± > 103.5 GeV and mgji > 100 GeV. The limit on m^^ is parametrized as a function 
of m^-^ as given by [36]; in case of a stau LSP, we take ruf^ > 94 GeV. For the light 
scalar Higgs, we apply the limits from LEP for the rn^""^ scenario given in [37] , taking into 
account a ~ 2 GeV theoretical error [38] . |^ 

We also take into account additional constraints from S-physics. For the branch- 
ing ratio of inclusive radiative B decay, we use the experimental result BR(6 — > 57) = 
(3.52 ± 0.23 ± 0.09) x 10"^ from HFAG [39], together with the SM theoretical prediction 
of BR(6 — > 37)^^^ = (3.15 ± 0.23) x 10""^ of [40]. Combining experimental and theoretical 
errors in quadrature, we require 2.85 < BR(6 57) x 10^ < 4.19 at 2cr. Another important 



^Moreover, the limits from direct squark and gluino searches at the Tevatron for rriq ~ rrig apply, in 
particular rrig > 392 GeV, but these are automatically fulfilled here. 
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constraint comes from the Bg decay into a pair of muons. Here we apply the 95% CL upper 
hmit BR(i?s f^~^f^') < 5.8 x 10~® from CDF [41]. Regarding the anomalous magnetic 
moment of the muon, we do not impose any limits but simply note that {g — 2)^ favours 
/i > 0. 

Last but not least, if the lightest neutralino is the LSP, we compare its relic density to 
the results from the 5-year WMAP data on the dark matter relic density, Vth"^ = 0.1099 ± 
0.0062 [42], although we do not impose this as a strict constraint. The values of BR(6 — > sj), 
BR{Bs /i"*"/^") and flh"^ are computed using the micrOMEGAs2 . 2 package [43]. 



6 Results for simplified boundary conditions 

Here we perform a first exploration of the parameter space using simplified boundary 
conditions in the matter sector according to [18]. More precisely, we assume that not only 
the first two generations but also the third-generation leptons and r.h. bottom are brane- 
localized. The top and l.h. bottom have a flat profile in the fifth dimension. The relevant 
kinetic functions then ar^ 

>Q3^n/3^|(T + T), (34) 



which leads to 



and 



2R 



2R1 + C' 

The setup for the first parameter scan is therefore as follows: 



(35) 
(36) 



We vary M1/2 from 100 and 1000 GeV and tan/? from 2 and 20. (For higher values 
of tan/5, the bottom Yukawa coupling would be no longer negligible.) 

We set mfj,^ = mg^ = My^{l + c')^ and At = M1/2; this requires the additional 
iteration on c' as detailed in Section [H All other sfermion soft terms are assumed to 
be zero at the GUT scale. 

We allow for all four sign combinations of sign(yu) = ±1 and en = ±1. 

/i and Bfi are determined from Eq. (29) at the EWSB scale, while m|^^ and m^^^ are 



determined from Eq. (31) at Mgut- 



For each point that gives correct EWSB, we check the mass limits from LEP as well 
as the constraints from BR(6 57) and BR(i?s At^/^~) given in Section [s] 



Figure |3] shows the result of this scan in the tan/5 versus M1/2 plane. As expected, correct 
EWSB is obtained only for two of the four possible combinations of sign(yu) and e^. In 
particular, it turns out that the two signs need to be equal. This is a consequence of 



^In this approximation yt — 54 is the only non- vanishing Yukawa coupling. Wc will however not enforce 
this in the numerical analysis. 
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Figure 3: Parameter points giving correct EWSB from a scan over M1/2 and tan/3 with 
simplified boundary conditions. The red, green and blue points have a neutralino, stau and 
selectron LSP, respectively. Small crosses denote points excluded by LEP, while open circles 
denote points excluded by B-physics constraints. The big full dots pass these constraints. 



the relation At = M1/2, in accord with the discussion in Section |4| Phenomenological 
bounds further constrain the parameter space. Points marked as small crosses in Fig. |3] are 
excluded by the mass bounds from LEP, while points shown as open circles are excluded 
by BR(6 57); the constraint from BR(i?s — >■ 11^ has no effect. The remaining big full 
points are phenomenologically viable. The different colours denote the nature of the LSP: 
red for a neutralino, blue for a selectron and green for a stau LSP. As one can see, most of 
the parameter space features a neutralino LSP, which is interesting in point of view of dark 
matter]^ As anticipated in Section [ij |/i| turns out to be large throughout the parameter 
space. Numerically we find |yu| ~ (2.5 — 3.5)Mi/2 for /i > and ~ (2.5 — 4)Mi/2 for 
/X < 0; in both cases the values at the high end are obtained for larger tan [3. The Xi is 
hence almost a pure bino, and the X2 ^'^d Xi almost pure winos. 

The projections onto the space of fundamental model parameters F^/2i?, F*^ and 
d are shown in Fig. |4j We observe that for both, /i < and /z > 0, there is a strong 
correlation between F"^ and F'^ /2R, with roughly F'^ ~ 3 x F^/2R. This comes from 
setting At = M1/2, which enforces e^j = sign(/i). It translates into a large value of F'^, 



because = e^fi + F'^ /2R from Eq. (pi). Nevertheless F'p is small enough so that 



contributions from anomaly mediation, being (9(F'^/87r^), can safely be neglected. 

It is particularly interesting to note that we find no valid spectra for which c' = 0. This 
also holds when considering points excluded by LEP constraints. In this sense our analysis 



"'^^Selectrons and smuons are taken to be mass-degenerate. Here and in the following we only refer to 
selectrons for simplicity, implicitly meaning "selectrons and smuons" . 

^^Alternative dark matter candidates would be gravitino or axino. A rough estimate for a no-scale radion 
Kahler potential K = -31og(T + T) gives mg/a > |F^/2i?|, while m^o ~ Q A Ml/2- In this case a gravitino 
LSP is only possible for c' < —0.6, which does not occur in our analysis. Moreover, we expect other 
contributions from hidden sectors to further increase 7713/2- An axino LSP is a valid option, but leads to 
a very different phenomenology, beyond the scope of this paper. 



14 




1000 1500 2000 2500 3000 

F'^nR [GeVl 



/i < 0, = -1 



2 




1500 2000 2500 3000 




1000 1500 

F^I2R [GeVl 



2000 




F'^/2R [GeVl 



1000 1500 

F'^/2R [GeVl 



2000 



Figure 4: Scatter plot of points which give a vahd spectrum solution in the F^/2i? vs. 
F"^ plane (top row) and in the F^/2i? vs. c' plane (bottom row). The red, green and blue 
points have a neutralino, stau and selectron LSP, respectively. Open circles denote points 
excluded by B-physics constraints. Points excluded by LEP are not shown. 
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Figure 5: Mass spectrum in the neutralino LSP region, passing LEP and B-physics con- 
straints, as a function of M1/2. The colour convention is as follows: red: Xi; green: fi, blue: 
Cij, dark blue: Cl, , dark green: f2, dark red: X2- 

confirms the result of [18], who did not include the effects of a Chern-Simons term and 
consequently did not find any viable parameter regions, except for extremely unnatural 
values for the gaugino masses far above our scan limits. At the same time it is important 
that c', which is an 0{1) parameter, never becomes large. 

Implications for collider phenomenology can be deduced from Fig. |5} which shows the 
neutralino and slepton mass spectrum in the neutralino LSP region. We see that the second- 
lightest neutralino X2 lighter chargino xf^ which are mainly winos {'m^± ~ mx^), 
are always heavier than en^i and fi (with the exception of a few points at yU < which have 
mg^ > m^o > rngj^). Note the clear separation of the selectron masses with mg^ < mg^ for 
H < 0, while for > we have mg^ ~ ^sr- The squark and gluino masses are not shown, 
but they are roughly rrig ~ rrig ~ (2 — 3)Mi/2- At the LHC, squarks and gluinos will hence 
be produced both as qq or gg pairs, and in qg associated production. Their decays are 
g — > qq^L, q^ — > qXi^ Ql ~^ I'xf or 1X2^ as in the mSUGRA scheme with large [44]. 
Moreover, the decays X2 ~^ ^^^l ~^ ^^^~Xi and X2 ~^ t^t^" t^t^Xi are always open 
and together have about 50% branching ratio; the other 50% go into neutrinos. This leads 
to the gold-plated same-flavour opposite-sign (SFOS) dilepton signature at the LHC [45], 
which allows to reconstruct sparticle masses. Furthermore, the decay of the lighter chargino 
always leads to a charged lepton, xf ~^ {^^^i or vd^) t^viXii giving rise to a large 
number of events with jets plus 1 hard lepton plus missing transverse energy, E^^^^. If 
combined with X2 ~^ ■■■ ~^ ^~^^~XiJ this leads to the rather clean trilepton signature (plus 
jets plus E^'^"). 

The scenario becomes even more predictive if we require that the neutralino LSP have a 
relic density in agreement with cosmological observations (assuming standard cosmology). 
Imposing the 3a upper bound from WMAP5, flh"^ < 0.1285, constrains M1/2 ^ 390 GeV 
with tan/3 > 11 for /i > 0. For /i < 0, it gives an upper limit on M1/2 which increases 
with tan/3, from M1/2 < 312 GeV at tan/3 = 12 to M1/2 < 920 GeV at tan/3 = 20. The 
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reason is that the LSP is almost a pure bino and has a small pair-annihilation cross section 
(s-channel Higgs exchange is not efficient in this scenario); in order to have a small enough 
relic density, the LSP needs to co-annihilate with another sparticle which is close in mass, 
typically the next-to-lightest SUSY particle (NLSP). This constrains the scenario to the 
region of small NLSP-LSP mass differences near the boundary to the slepton LSP region, 
which is reahzed for /i < up to large M1/2 (depending on tan/3), but for /i > only 
at small M1/2, cf. Figs. |3] and |5j Note, however, that this is a direct consequence of the 
simplified assumptions for the matter sector. 



7 Realistic sfermion soft terms 

In this section we explain how improved sfermion soft terms can be obtained if we model 
the matter sector as in the Burdman-Nomura model [12]. The third generation matter 
fields arise from the mixing of brane and bulk fields. Bulk fields with flat proffie have 
Yukawa couplings determined by the 5d gauge coupling. Non-trivial bulk profiles cause a 
reduced overlap with the Higgs wave function and hence smaller Yukawa couplings. Thus, 
using both bulk masses and mixing angles we can obtain realistic values for yt, yi, and yr- 

For the third generation quarks in particular, we introduce a 5d bulk hypermultiplet 
{W, W^} in the 20 of SU(6) containing as 4d zero modes the right-handed top quark super- 
field and a weak doublet, and another bulk hypermultiplet {"D, P^} in the 15 containing 
the right-handed bottom quark and a second doublet. We give these fields bulk masses 
Mu and Ma- Furthermore, brane-localized superfields must be introduced to decouple un- 
wanted massless fields. They couple to the doublet components of both the U and V fields, 
leaving a single massless quark doublet instead of the two we were starting with. This 
effect is parametrized by a mixing angle 0q. 

Similarly, leptons descend from two 5d bulk hypermultiplets, {S,S'^} in the 15 and 
{N',N''^} in the 6. In analogy with the quark sector this leads to three more model pa- 
rameters, two bulk masses Me and Mn and a mixing angle 0^. For details of the model, 
in particular for the proper choice of boundary conditions, brane fields and bulk-brane 
couplings, we refer to [12]. 

The kinetic functions are computed by integrating the zero-mode proffies over the fifth 
dimension, replacing its radius by (T + T)/2. This gives 

>^^3 = ;^fl-e-'^™M, (37) 



Yq, = ^ (1 - e-(^+^)l^^^l) sin2(0g) + ^ (1 - e-(^+^)l^^^l) cos\<Pq), (38) 

r.3^^(l-e-™). (39) 

The kinetic functions for the lepton sector are obtained in the same manner, and are 
given by the same expressions with the obvious parameter replacements. The soft masses 



and A-terms are then derived from Eqs. (25) - (27). We refrain from giving closed-form 



expressions for them, since these are rather cumbersome and not very illuminating. 

The parameters Mu, Ma, Mn, Me, (pq and 0l cannot be chosen entirely freely, because 
they also have to account for the proper physical values of the Yukawa and gauge couplings. 
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Since the Higgs wave function normahzation is just given by (Yh) = Vs'I ^^"^ i'^ particular 
is independent of c', the relations given in [12] apply |^ 

T^R\Md\ 



yt = sin(0Q) I 94, 



sm 



sinh ttR\Mu\ 

7lR\Mn\ 

sinh ttR\M„ 



Vb 



94, Vt 



cos 



cos 



sinh 'KR\Md 

TlR\Me\ 

sinhvri^lM^ 



94 



94 



(40) 
(41) 



In the numerical analysis, in order to avoid additional model dependence from the 
unknown neutrino sector, we will assume that M„ is large enough not to contribute to the 
stau soft terms. We also introduce a Majorana mass term for the right-handed neutrinos on 
the y = brane as in [12]. Since M„ is large, the neutrino wave function will be strongly 
localized towards the y = nR brane, resulting in an exponentially suppressed Yukawa 
coupling and a doubly exponentially suppressed Majorana mass. The suppression factors 
will cancel out in the see-saw formula for the lighter neutrino mass eigenstate, leading to 
the same lighter neutrino mass as in the standard see-saw mechanism. The heavier neutrino 
mass, on the other hand, will be lowered by a factor ~ g-47ri?|M„| ^j^]-^ respect to the GUT 
scale. This may be beneficial for leptogenesis [46]. 



It is instructive to see how Eq. (40) constrains the possible ranges of squark soft terms. 



In the remainder of this section we will therefore give some estimates of the bounds on the 
squark masses and trilinear couplings. 

For tan /3 ~ 5 — 50, the relevant GUT-scale Yukawa couplings take values 0.5 < < 0.6 
and 0.02 ^ ?/& < 0.3. We also know that the gauge couplings unify at ~ 0.7. 

To reproduce the top Yukawa coupling, we must have tan^g > 1 by Eq. (40). The 
small ratio yb/yt can then be generated either by choosing tan^g to be large, or choosing 
\Mii\ > |M„|, or by a combination of these. The relation between M^, M^, 0q and the 
Yukawa couplings is illustrated in Fig. |6} We note that for given yt and yt, the allowed 
range for the mixing angle 0q is 



[ aicsin^yt/g^) , arccos(?/;,/ g^)]. 



(42) 



For estimating the size of the squark-mass parameters, let us consider two limiting cases: 

• If the difference between yt and yb is mainly due to the different bulk masses, then 
tan0Q 1. This corresponds to the far left region of Fig. [6} In that case sin$Q ^ 
l/-\/2 already accounts for the ratio 7/4/(74 ~ 0.7in Eq. (40). The top Yukawa coupling 
should thus not receive much additional suppression from large bulk masses, hence 
we need |M„| <C 1/i?. Expanding Eq. (37) and retaining only the leading term, we 



reproduce Yu^ as in Eq. (34) 



Yr 



Us 



(T + T) 



m 



Us 



2R 



(43) 



On the other hand, R\Md\ must be sizeable to obtain an appropriately suppressed 
yb, cf. Fig. |6j With Eq. (39), m|,^ turns out to be 

2R 



m 



'nR\Md\ 
sinh(7ri?|Md|) 



4y, 



(44) 



^^Note that our conventions for 0q and slightly differ from those of [12]. 
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Figure 6: Values of -R|M„| and R\Mfi\ as function of the mixing angle 0q for various values 
of yt (full red lines) and i/b (dashed blue lines). Note that yt gives the lower and yi, the 
upper bound of the allowed range of (pQ. 



Finally, the quark doublet soft mass-squared uiq^ obtained from Eq. (38) is numeri- 
cally 



(0.7-0.8) X 



2R 



(45) 



If tan0Q ^ 1, i.e. sin^g ^ 1 (which is the case in the far right region of Fig. [6]), 
then the ratio yt/g^ ~ 0.7 is mainly due to a sizeable bulk mass M„. Numerically, we 
need R\Mu\ ~ 0.3 — 0.5. Therefore we should use the full expression for Yjj^, rather 
than just the leading term: 



m 



Us 



7lR\Mu\ 



sinh(7ri?|M„|) 



2R 



(0.5-0.8) X 



2R 



(46) 



Dropping the cos^ 0q piece in Yq^ and setting sin ( 
pression for Yq^ and eventually rriQ^: 



1, we obtain the same ex- 



m 



nR\M^\ 



sinh(7ri?|M„|) 



2R 



(0.5-0.8) X 



2R 



(47) 



As is evident from Fig. |6| if yi, is to remain finite, tan 0q cannot become arbitrar- 
ily large. In any case, this limit requires very small yt,. The constraints on \Mii\R 
are rather weak, although smaller \Md\R is somewhat favoured in order not to get 
additional yt suppression. Hence 



2 

D3 ~ 



2R 



(4J 
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In the end we expect the squark masses-squared to he somewhere in between these two 
extremes: 



0.5 X 



2R 



< 



< 



2R 



< < 



2R 



In order to obtain limits on At, we can make the same case distinction: 



(49) 



for tan i 



1 and small |Mu|, we get 




1 + 2c' 



1 + 



2TiR\Md\ (1 + e 



-27r/?|MdP 



2ti R\Md\ + 1 - e-2^«l^^'*l 



(50) 



for sini 



1, we obtain 



At 



2R 



27ri?|M, 



1 + 2c' 



1 + c' exp(27ri?|M„|) - 1 



(51) 



Numerically, 



At 



2R 



1 + 2c' 
1 + c' 



a 



(52) 



where 0.3 < a < 2, with a = 2 corresponding to the first of the above two cases (with 
R\Md\ ~ 1), and a = 0.3 to the second (with R\Mu\ = 0.5). Evidently At can take 
a wide range of values, significantly departing from the simplified case of Section [6] In 
particular it can become large and negative, which will be of relevance in the next Section. 
A similar statement turns out to be true for Af,, for which we find an analogous estimate 
with < a < 1.4. 



8 Results for realistic sfermion soft terms 

Let us finally investigate to what extent the phenomenological features found in Section [6] 
remain valid when invoking realistic stop, sbottom and stau parameters derived from the 
Burdman-Nomura model. The six new parameters Mu, Ma, Mn, Me, 0q, (pL are subject 



to four constraints, since they are related to the Yukawa couplings according to Eqs. (40) 



and (41). As detailed above we assume that M„ is large enough not to affect the stau soft 
terms. This corresponds to a neghgible neutrino Yukawa coupling, and we do not need to 
worry about lepton flavour violation [47]. The precise value of M„ is irrelevant. (If M„ did 
contribute to the stau soft terms, its main effect would be to increase m^g, thus rendering 
the staus heavier, but leaving the overall picture intact.) We are therefore left with five 
parameters, M^ ^ e and 0q,l, and three constraints from yt, yb and y^. 

We choose 0q and 0^ as the two independent new parameters and scan the parameter 
space as in Section |6| with the following modifications: 

• We vary M1/2 from 100 to 1000 GeV, 0q from 7r/4 to 7r/2, and 0l from to 7r/2. 
For tan/5, we consider three distinct values, tan/5 = 10, 20, and 30, in order to avoid 
excessive computing times. 
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For each point, the bulk masses Mude are computed from the GUT-scale gauge and 



Yukawa couphngs 5(4, i/t, i/b, Ut by numerically inverting Eqs. (40) and (41). They 



then serve as input in the kinetic functions Eqs. (37)-(39), and the analogous expres 



sions for the leptons, from which the sfermion soft masses and A-terms are obtained 



according to Eqs. (25)-(27). The soft terms of the first and second generation are 



again assumed to be zero at the GUT scale. 

The result of this scan is shown in Figs. [7] and [8] for the two signs of /i. For better 
readability, we only show M1/2 in steps of 200 GeV, although the scan had a much finer 
grid. Contrary to the case of simplified boundary conditions, now and e^j need to be of 
opposite sign. The reason is that now At turns out to be negative at the GUT scale (cf. the 
discussion in Section lil). 



It is interesting to see how the mixing angles 0q and (pL influence the nature of the 
LSP. (pL determines the size of the stau parameters. Since it is constrained by the tau 
Yukawa coupling it can only vary over a sizable range if tan (3 is large. The reason is that 
Ar is generically large, leading to a charge-breaking minimum if uil^^Es are too small. 
Thus for tan/3 ~ 10, (pi is close to 7r/2 and the staus are rather heavy compared to the 
selectrons. For larger tan/5 (i.e. larger (p^ can be small and the fi can become the 
LSP, corresponding to the green points in Figs. [7] and [8| Note, however, that the stau LSP 
region is highly constrained by direct mass bounds and -B-physics, and that a stable LSP 
is excluded by cosmology. 0^ also has some effect on the selectron masses through RG 
evolution, but this is much less pronounced. 

The angle 0q, on the other hand, determines the size of the stop and sbottom pa- 
rameters. Through RG evolution it also influences the slepton masses, in particular mg^: 
larger 0q leads to a larger m^g, which in turn decreases mg^. In Figs. [7] and [8] one can see 
clearly that for increasing 0q, the cr eventually becomes the LSP. This behaviour can be 
understood easily from the U(1)y -D-term contribution to the evolution of the scalar soft 
masses mf [48]. At one loop 

d 2 6 glYi 
where is the weak hypercharge and 

^ = i^H2 - ^Hi) + Tr {ml - 2ml + ^D + mR- ^1) (54) 

with the trace running over generations. Since 5* is an RG invariant, it simply causes a 
shift of the low-scale masses by Amf ^ —(0.052) Squt [49] with respect to the values 
they would have had for 5 = 0. Here Sgut is the value of 5* at Mqut- For simplified 
boundary conditions, we had Squt = ~m1^. With Y^.^ = 1 and Y^^ = —1/2, making Squt 
less negative obviously lowers mg^ and increases mg^ (note also that the effect for the 
left-chiral state is only half the size of that for the right-chiral one). Moreover, comparing 
mg^ ~ (0.39 Mi/2)^ — 0.052 S'gut to m^o ^ 0.43 M1/2, we understand why the cr eventually 
becomes the LSP. 



The projections onto the underlying model parameters F'^ /2R, F"^ and c' are shown 
in Fig. [9} Since here we need sign(/i) = —ej^ to obtain a valid spectrum, F*^ now turns out 
to be small and can even be zero. Contributions to the soft terms from anomaly mediation 
are therefore completely negligible. Moreover, we find a somewhat smaller range for the c' 
parameter, roughly 0.5 < c' < 1.2, as compared to 0.5 < c' < 3 for simplified boundary 
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conditions. The important point, however, is that c' remains non-zero. We conclude that 
the Chern-Simons term is indeed essential to achieve correct EWSB. 



Let us now turn to the implications for collider phenomenology. We again focus on 
the neutralino LSP region. The mass spectrum in this region, taking into account the 



constraints from LEP and from B-physics, is depicted in Fig. [TOj As one can see, there 
is a definite mass ordering m-± ~ m^o > TJip, > mp„ > m^o. The fo turns out to be 
heavier than the while the fi can be lighter than the X2, and for small 0^ also lighter 
than the selectrons, cf. the above discussion of the mixing-angle dependence. This gives a 
picture that is qualitatively similar to the simplified case discussed in Section |6] the main 
difference lies in the masses and mass ratios of the sleptons. For the squarks, this effect 
of non-universality — on the one hand the splitting of the third generation from the first 
and second generations due to non-zero ^^3, {73,1)3, on the other hand the splitting of left- 
and right-chiral states due to non-zero S — is much less pronounced, because the running 
of the squark mass parameters is mainly driven by M3. The squark and gluino masses are 
hence again about rrig ^ rrig ^ (1.7 — 2.5)Mi/2- The masses of the higgsino-like neutralinos 
and chargino are given by and lie above rrig. 

It is also remarkable that now the neutralino relic density can vary over a large range, 
because of the extra parameters 0q and 0^. This is illustrated in Fig. 11 for the example 
of Ml/2 = 500 GeV and two values of tan/3 for each sign of /i. For /x > 0, we take 
tan P = 10 and 30; for /i < 0, we take tan /3 = 10 and 20 since higher values are too tightly 
constrained. The figure compares the neutralino relic density Qh"^, as a function of 0q 
and (pL, with the WMAP5 observation at 3a. In the orange regions flh"^ is too low, which 
would require other constituents of dark matter in addition to the neutralino. In the brown 
regions, on the other hand, Qh"^ is too high (at least within standard cosmology; it could 
be viable if there was, e.g., additional entropy production after freeze-out). The minimal 
and maximal values found are flh^ ~ 6 x 10~^ and 0.9, respectively. In the red band in 
between, however, 0.0913 < Qh"^ < 0.1285 agrees within 3a with the value measured by 
WMAP5. The reason is that here the mass difference between the LSP and NLSP (or 
co-NLSPs) is just right to make co-annihilation processes efficient enough, but not too 
efficient, to obtain Qh"^ ~ 0.1. To be precise, in the red bands of Fig. 11 we typically have 
Am = rrieR — m^o ^7 — 10 GeV. An exception is yU > 0, tan/3 = 30 and small (pL, where 
the fi becomes light and also contributes to co-annihilations, such that Am ^ 20 GeV 
is needed; this leads to the red band bending down towards lower (pq. Sample spectra of 
five representative points, indicated as points A-E in Fig. [TTJ are given in Tables [T] and [2j 
Table [T] gives GUT and EWSB scale parameters, and Table [2] lists the resulting masses 
together with B-physics observables, the neutralino relic density and the neutralino-proton 
scattering cross section for direct detection. The possibility to tune the NLSP-LSP mass 
difference by adjusting 0q and 0l and to obtain the correct relic density persists also for 
other values of Mi/2- 

To summarize, the expected LHC phenomenology is as follows: 



• Squarks and gluinos with masses up to about 2 TeV will be abundantly produced 
at the LHC, both as qq or gg pairs, and in qg associated production. They decay as 
~9 ^ Q'lR,L^ Qr ^ QX^ (~ 100%), q^ ^ q'xf (~ 65%) or qxl (~ 30%). 

• The decay X2 — ^ e^e^ — > e^e~Xi is always open and has a sizable branching ratio 
(~45% for points A-D, 35% for point E). This leads to a rather large rate for the gold- 
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Figure 7: Points which lead to correct EWSB from a scan over M1/2, tan/3, and 
for /X < 0, = +1 and sfermion soft terms determined according to the Burdman- 
Nomura model Small crosses denote points excluded by LEP, while open circles denote 
points excluded by 5-physics constraints. Points passing these constraints are shown as 
big full dots. The colours denote the nature of the LSP: red, green and blue points have a 
neutralino, stau and selectron LSP, respectively. 
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Figure 10: Mass spectrum in the neutralino LSP region, passing LEP and B-physics con- 
straints, as a function of M1/2. From bottom to top: (red), tn (blue), (dark blue), fi 
(green), X2 (dark red), f2 (dark green) and g (purple). 



plated SFOS dilepton signature. In parts of the parameter space, also X2 ~^ t^t^ — > 
r~r+x5 can be kinematically allowed; c.f. point E, where it has 22% branching ratio. 
Decays into Z, /i, or are negligible because the X2 is almost a pure wino. 

• The decay of the Xi" always leads to a charged lepton, — » t^Vijvil^ — >• i^uiXi, 
giving rise to a large number of events with jets plus one hard lepton plus EJp^^^. 
If combined with X2 ~^ ■■■ ~^ ^^^~Xi on the other side of the event, it leads to the 
rather clean trilepton signature (plus jets plus E^^^^). 

• The higgsino states X34 and xt have masses around or above the gluino mass and 
are hence too heavy to be studied at the LHC. 



Overall the scenario resembles the mSUGRA/CMSSM case with small mo, or the case 
of Higgs boson exempt no-scale supersymmetry (HENS) [49]. An important difference are 
the sizeable third-generation high-scale soft terms which our construction predicts. Ways 
to distinguish between the different models include, e.g., the rate of leptonic events, which 
is expected to be higher in our scenario as compared to the mSUGRA case with the same 
Ml/2. Other distinctive features are the ratios of left- and right-chiral slepton masses and 



the non-universality of the third generation,^ Note, however, that the cr does not couple 
to the wino-like X2 xf and hence does not appear in decay chains at the LHC. The 
cr is therefore best studied in e~^e~ collisions, as are the staus if they are too heavy to be 
produced in X2 decays. 

Last but not least, a decisive test of GHU requires the precise measurement of the 
complete spectrum, including stops, sbottoms, heavy Higgs bosons and higgsinos, such 
that the SUSY Lagrangian parameters can be extracted and a bottom-up evolution along 



^^We leave a detailed study of characteristic mass ratios and model footprints for future work. 



26 




Figure 11: Dependence of the neutralino relic density on the mixing angles 0q and 0l, for 
Ml/2 = 500 GeV and various values of tan/?. In the red bands, lies within Sex of the 
WMAP5 observation, 0.0913 < VtK^ < 0.1285. In the orange regions, Vth^ < 0.0913 is too 
low, while in the brown regions Qh"^ > 0.1285 is too high. Also indicated are the sample 
points A-E. 



the lines of [50] performed. This can only be achieved at a (multi-)TeV e~^e linear collider 
with a very good beam performance. 



9 Conclusions 



We have investigated SUSY grand unified models with gauge- Higgs unification (GHU). 
A particularly interesting class of such models are 5d orbifold GUTs and heterotic string 
models which admit a 5d orbifold GUT limit. With the natural assumption of radion 
mediation, GHU models are quite predictive as far as the Higgs sector is concerned. The 
GUT-scale Higgs mass parameters are subject to the GHU relations, and are also tied 
to the gaugino mass. Despite these strong constraints, models of this type can be fully 
realistic, as we have shown. If the effects of a Chern-Simons term (which is generically 
present in 5d models) are taken into account, one finds regions in the parameter space 
which lead to proper electroweak symmetry breaking and satisfy the experimental bounds 
from direct Higgs and superpartner searches, rare decays and cosmology. We demonstrated 
this by using a variation of a 5d SU(6) orbifold GUT model due to Burdman and Nomura 
as a concrete example. We gave detailed expressions for the soft SUSY breaking parameters 
in terms of the fundamental model data, including the Chern-Simons term. 

Using the high-scale relations between soft terms and estimates of running effects, 
we discussed qualitatively which parts of the parameter space might be promising. We 
then presented a detailed numerical analysis of the corresponding RGEs. This analysis was 
done in two parts, the first for a simplified model of the sfermion sector, and the second 
treating the relevant sfermion contributions properly as in the Burdman-Nomura model. 
The latter part of the analysis, while more realistic, is more involved because it depends 
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on more parameters. In both cases we indeed find viable solutions to the RGEs, satisfying 
all present experimental constraints. A non-zero Chern-Simons term is essential to get a 
valid spectrum. 

We extracted some characteristic experimental signatures of this class of models, which 
will be tested at the LHC In particular, selectrons are generically predicted to be lighter 
than the leading to a rather large rate for same-flavour opposite-sign dileptons over the 
whole parameter space. Higgsinos, on the other hand, are expected to be heavy, presumably 
beyond the reach of the LHC. Characteristic mass ratios could be tested in detail at a future 
e'^e~ linear collider. 

The LHC will have the potential to narrow down the allowed region in the parameter 
space of GHU models significantly, or to rule them out. This applies even more to a future 
linear collider. It would be worthwhile to study in detail how well the scenario discussed 
here could be reconstructed at the LHC and a linear e~^e~ collider, thereby testing the 
GHU relation. To this end we proposed a set of benchmark points which may be useful for 
Monte Carlo simulations. 
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Point 


A 


B 


C 


D 


E 


Ml/2 


500 


500 


500 


500 


500 


tan/5 


10 


20 


10 


30 


30 


sign(/i) 


-1 


-1 


+ 1 


+ 1 


+ 1 


en 


+ 1 


+ 1 


-1 


-1 


-1 


<pQ 


1.3011 


1.1503 


1.3486 


1.1582 


1.1027 


<Pl 


1.2376 


1.0314 


1.3329 


0.8889 


0.1437 


c' 


0.5712 


0.6118 


0.7577 


0.6686 


0.6811 


F^/2R 


785.6 


805.9 


878.9 


834.3 


840.9 




-107.7 


-120.6 


-38.9 


-114.6 


-109.9 


Parameters at Mqut 




-1178.9 


-1232.4 


1296.6 


1283.2 


1291.0 


B 


-344.0 


-381.3 


-298.0 


-393.2 


-390.2 




564.8 


615.2 


622.9 


634.9 


657.3 




239.5 


371.2 


279.5 


413.9 


369.2 


mq^ 


555.5 


592.4 


615.9 


613.5 


626.3 


At 


-812.6 


-793.7 


-980.0 


-842.2 


-823.1 


A, 


-932.3 


-923.6 


-1105.6 


-966.3 


-978.8 




226.7 


304.9 


364.1 


432.6 


278.8 




161.8 


269.4 


225.2 


407.9 


278.5 


Ar 


-1055.8 


-1073.7 


-1216.8 


-1070.4 


-1149.6 


Parameters at Mewsb 




-1217.6 


-1238.8 


1342.0 


1275.3 


1282.9 


B 


-46.8 


-16.7 


45.9 


11.4 


10.6 


2 


554195 


405781 


598972 


428198 


401063 


ml 


2022 


-2969 


2560 


-3429 


-3467 



Table 1: Parameters of sample points A-E in Fig. 11 Dimensionful quantities are in GeV, 
B = Bfi/fi. 
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Point 


A 


B 


C 


D 


E 


Al 


210.1 


210.3 


208.1 


208.7 


208.7 




389.3 


389.5 


399.2 


400.4 


400.3 




1219.7 


1240.7 


1332.0 


1265.2 


1272.8 


A4 


1220.3 


1241.9 


1335.1 


1267.7 


1275.3 


m-± 

Xi 


389.3 


389.4 


399.2 


400.4 


400.3 


-± 

A 2 


1222.8 


1211.1 


1335.1 


1268. 1 


1276.0 






'^98 1 
OZo. i 


'^97 


OZO.O 


OZo. i 




Zio.O 


91 7 n 


91 


917 1 


998 
zzo.o 




9Qp; Q 

zyo.y 


^^99 n 




^^87 7 


99c; A 
ZZoA 


"''T2 


372.6 


441.2 


438.0 


549.5 


457.0 


TTI.r, 


318.4 


318.3 


318.1 


317.3 


313.4 




354.3 


408.0 


386.2 


499.0 


398.4 












1 n/i 1 Q 


TYl ~ 


1003.5 


1003.2 


1000.0 


999.6 


997.6 


771 J 


1049.0 


1048.7 


1045.1 


1044.8 


1044.9 


777 r 


1005.9 


1005.6 


1001.9 


1001.8 


1002.3 


777 r 


948 


971 1 


955 7 


971 4 


983 5 


777r 


1147.0 


1155.0 


1187.9 


1167.6 


1175.2 




1022.3 


1016.7 


1029.9 


1021.7 


1007.9 




1108.4 


1119.6 


1137.6 


1130.6 


1135.5 


mg 


1155.5 


1156.5 


1154.2 


1154.7 


1155.0 




115.0 


116.4 


117.2 


117.3 


116.8 




762.4 


658.5 


770.7 


637.4 


635.9 


rriA 


761.6 


658.5 


770.8 


637.6 


635.9 


mH± 


766.7 


663.9 


775.1 


642.9 


641.0 


BR{b S7) 


3.70 X 10-4 


4.16 X 10-4 


3.20 X 10-4 


2.89 X 10-4 


2.91 X 10-4 




2.89 X 10-9 


2.10 X 10-9 


3.06 X 10-9 


6.76 X 10-9 


6.68 X 10-9 




0.110 


0.108 


0.110 


0.108 


0.106 


(t(xp)^^ [pb] 


2.92 X 10-11 


1.39 X 10-1° 


1.01 X 10-1° 


2.90 X 10-1° 


2.89 X 10-1° 



Table 2: Masses (in GeV), B-physics observables, relic density and spin-independent 
neutralino-proton scattering cross section for points A-E. 



30 



References 



[1] H. Georgi and S. L. Glashow, "Unity Of All Elementary Particle Forces," Phys. Rev. 
Lett. 32 (1974) 438; 

H. Georgi, H. R. Quinn and S. Weinberg, "Hierarchy Of Interactions In Unified Gauge 
Theories," Phys. Rev. Lett. 33 (1974) 451. 

[2] S. Dimopoulos, S. Raby and F. Wilczek, "Supersymmetry And The Scale Of Unifica- 
tion," Phys. Rev. D 24 (1981) 1681. 

[3] D. B. Fairlie, "Higgs' Fields And The Determination Of The Weinberg Angle," Phys. 
Lett. B 82 (1979) 97; 

N. S. Manton, "A New Six-Dimensional Approach To The Weinberg-Salam Model," 
Nucl. Phys. B 158 (1979) 141; 

Y. Hosotani, "Dynamical Mass Generation By Compact Extra Dimensions," Phys. 
Lett. B 126 (1983) 309. 

[4] W. Buchmiiller, K. Hamaguchi, O. Lebedev and M. Ratz, "Supersymmetric stan- 
dard model from the heterotic string," Phys. Rev. Lett. 96 (2006) 121602 [arXiv:hep^ 



ph/0511035 l and Nucl. Phys. B 785 (2007) 149 |arXiv:hep-th706 06187l; 
O. Lebedev, H. P. Nilles, S. Raby, S. Ramos- Sanchez, M. Ratz, P. K. S. Vaudrevange 
and A. Wingerter, "A mini-landscape of exact MSSM spectra in heterotic orbifolds," 
Phys. Lett. B 645 (2007) 88 |a rXiv:hep-th/0611095] . 

[5] H. P. Nilles, S. Ramos- Sanchez, M. Ratz and P. K. S. Vaudrevange, "From strings to 
the MSSM," .arXiv:0806.3905| [hep-th]. 

[6] Y. Kawamura, "Triplet-doublet splitting, proton stability and extra dimension," Prog. 
Theor. Phys. 105 (2001) 999 |arXiv:hep-ph/0012125|; 

G. Altarelli and F. Feruglio, "SU(5) grand unification in extra dimensions and proton 
decay," Phys. Lett. B 511 (2001) 257 |arX iv:hep-ph/0102 301|; 

L. J. Hall and Y. Nomura, "Gauge unification in higher dimensions," Phys. Rev. D 
64 (2001) 055003 |arXiv:hep-ph/0103125l ; 

A. Hebecker and J. March-Russell, "A mi nimal S(1)/(Z(2) x Z'(2)) orbifold GUT," 
Nucl. Phys. B 613 (2001) 3 |arXiv:hep-ph/0106166|; 

T. Asaka, W. Buchmiiller and L. Covi, "Gauge unification in six dimensions," Phys. 
Lett. B 523 (2001) 199 |arXiv:hep-ph/ 0108021|; 

L. J. Hall, Y. Nomura, T. Okui and D. Tucker-Smith, "SO(IO) unified theories in six 
dimensions," Phys. Rev. D 65 (2002) 035008 [ arXiv:hep-ph/0108071| . 

[7] V. S. Kaplunovsky, "Mass Scales Of The String Unification," Phys. Rev. Lett. 55, 
1036 (1985). 

[8] L. E. Ibariez and D. Liist, "Duality anomaly cancellation, minimal string unification 
and the effective low-energy Lagrangian of 4-D strings," Nucl. Phys. B 382 (1992) 
305 (aBGv:hepT h/9202046 1 . 

[9] E. Witten, "Strong Couphng Expansion Of Calabi-Yau Compactification," Nucl. 
Phys. B 471 (1996) 135 |arXiv:hep-th/9602070j . 



31 



[10] A. Hebecker and M. Trapletti, "Gauge unification in highly anisotropic string com- 
pactifications," NucL Phys. B 713 (2005) 173 |a rXiv:hep-th/041113 1]. 

[11] B. Dundee, S. Raby and A. Wingerter, "Reconcihng Grand Unification with Strings 
by Anisotropic Compactifications," arXiv:0805.4186 [hep-th]; 

B. Dundee and S. Raby, "On the string couphng in a class of stringy orbifold GUTs," 
| arXiv:0808.0992 |! [hep-th]. 

[12] G. Burdman and Y. Nomura, "Unification of Higgs and gauge fields in five dimen- 
sions," Nucl. Phys. B 656 (2003) 3 |arXiv:he p-ph/0210 257). 

[13] Y. Nomura, D. Poland and B. Tweedie, "Holographic Grand Unification," JHEP 0612 
(2006) 002 | arXiv:hep-ph/0605014 |. 

[14] H. M. Lee, H. P. Nilles and M. Zucker, "Spontaneous localization of bulk fields: The 
six-dimensional case," Nucl. Phys. B 680 (2004) 177 | arXiv:hep-th/0309195] . 

[15] W. Buchmiiller, C. Liideling and J. Schmidt, "Local SU(5) Unification from the Het- 
erotic String," JHEP 0709 (2007) 113 |arXiv:0707.1651 [hep-ph]]; 
W. Buchmiiller, R. Catena and K. Schmidt-Hoberg, "Small Extra Dimensions from 
the Interplay of Gauge and Supersymmetry Breaking," Nucl. Phys. B 804 (2008) 70 
|arXiv:0803. 4501 [hep-ph]]. 

[16] P. Hosteins, R. Kappl, M. Ratz and K. Schmidt-Hoberg, "Gauge-top unification," 
arXiv:0905.3323 [hep-ph]. 

[17] N. Haba and Y. Shimizu, "Gauge-Higgs unification in the 5 dimensional E(6), E(7), 
and E(8) GUT s on orbifold," Phy s. Rev. D 67 (2003) 095001 [Erratum-ibid. D 69 
(2004) 059902] |arXiv : hep-ph702l2 166|; 

I. Gogoladze, Y. Mimura and S. Nandi, "Gauge Higgs unification on the left-right 
model," Phys. Lett. B 560 (2003) 204 [arXiv:hep-ph/03010l4|; 

I. Gogoladze, Y. Mimura and S. Nandi, "Unification of gauge, Higgs and matter in 
extra dimensions," Phys. Lett. B 562 (2003) 307 | arXiv:h ep-ph/0 302176 ]; 
A. Hebecker and M. Ratz, "Group-theoretical aspects of orbifold and conifold GUTs," 
Nucl. Phys. B 670 (2003) 3 iarXiv:hep-ph /0306049]; 

I. Gogoladze, T. Li, V. N. Senoguz and Q. Shafi, "Non-canonical MSSM, unifica- 
tion, and new particles at the LHC," Phys. Rev. D 74 (2006) 126006 [ar Xiv:hep^ 
pi/0608181] ; 

Y. Mimura and S. Nandi, "Three family unification in higher dimensional models," 
arXiv:0902.4275| [hep-ph] ; 

Y. Mimura, "Flavor Symmetry in Gauge-Higgs-Matter Unified Orbifold GUTs," 
[a"rXiv:0903.1875 [hep-ph]. 

[18] K. w. Choi, N. y. Haba, K. S. Jeong, K. i. Okumura, Y. Shimizu and M. Yamaguchi, 
"Electroweak symmetry breaking in supersymmetric gauge-Higgs unification models," 
JHEP 0402 (2004) 037 [arXiv:hep-ph/0312178] . 

[19] G. Lopes Cardoso, D. Liist and T. Mohaupt, "Moduli spaces and target space duality 
symmetries in (0,2) Z(N) orbifold theories with continuous Wilson lines," Nucl. Phys. 
B 432 (1994) 68 [arXiv:hep-th/9405002] . 



32 



[20] I. Antoniadis, E. Gava, K. S. Narain and T. R. Taylor, "Effective mu term in super- 
string theory," Nucl. Phys. B 432 (1994) 187 | arXiv:hep-th/9405 024|. 

[21] A. Brignole, L. E. Ibafiez, C. Mufioz and C. Scheich, "Some Issues In Soft Susy 
Breaking Terms From Dilaton / Moduli Sectors," Z. Phys. C 74 (1997) 157 [arXiv:hep^ 
|ph/9508258]; 

A. Brignole, L. E. Ibafiez and C. Mufioz, "Orbifold-induced mu term and electroweak 
symmetry breaking," Phys. Lett. B 387 (1996) 769 |arXiv:hep^h/9607405]; 
A. Brignole, L. E. Ibafiez and C. Mufioz, "Soft supersymmetry-breaking terms from 
supergravfty and superstring models," arXiv:hep-ph/9707209 . 

[22] T. Kobayashi and K. Yoshioka, "Kaluza-Klein mediated supersymmetry breaking," 
Phys. Rev. Lett. 85 (2000) 5527 |arXiv:hep-ph/0008069|; 

Z. Chacko and M. A. Luty, "Radion mediated supersymmetry breaking," JHEP 0105 
(2001) 067 |arXiv:hep-ph/0008103). 

[23] G. F. Giudice and A. Masiero, "A Natural Solution to the mu Problem in Supergravity 
Theories," Phys. Lett. B 206 (1988) 480. 

[24] V. S. Kaplunovsky and J. Louis, "Model independent analysis of soft terms in ef- 
fective supergravity and in string theory," Phys. Lett. B 306 (1993) 269 ar Xiv:hep^ 
[th/93030 40|. 

[25] A. Hebecker, J. March-Russell and R. Ziegler, "Inducing the /i and the Bfi Term by 
the Radion and the 5d Chern-Simons Term," |arXiv:0801.4101_ [hep-ph]. 

[26] N. Arkani-Hamed, A. G. Cohen and H. Georgi, "Anomalies on orbifolds," Phys. Lett. 
B 516 (2001) 395 |arXiv:hep-th/0103135]; 

C. A. Scrucca, M. Serone, L. Silvestrini and F. Zwirner, "Anomalies in orbifold field 
theories," Phys. Lett. B 525 (2002) 169 |arXi v:hep-th/0110073 |; 
R. Barbieri, R. Contino, P. Creminelli, R. Rattazzi and C. A. Scrucca, "Anomalies, 
Fayet-Iliopoulos terms and the consistency of orbifold field theories," Phys. Rev. D 
66 (2002) 024025 |a rXiv:hep-th /0203039i. 

[27] C. A. Scrucca and M. Serone, "Anomalies in field theories with extra dimensions," 
Int. J. Mod. Phys. A 19, 2579 (2004) | arXiv:hep-th/0403163| . 

[28] A. Djouadi, J. L. Kneur and G. Moultaka, "SuSpect: A Fortran code for the super- 
symmetric and Higgs particle spectrum in the MSSM," Comput. Phys. Commun. 176 
(2007) 426 |arXiv:hep-ph/0211331], 

|http : //www ■ Ipta. univ-montp2 . f r/users/kneur/Suspect/ 

[29] D. Martf and A. Pomarol, "Supersymmetric theories with compact extra dimensions 
in N = 1 superfields," Phys. Rev. D 64 (2001) 105025 |arXiv:hep-th/010 6256]. 

[30] N. Arkani-Hamed, T. Gregoire and J. G. Wacker, "Higher dimensional supersymmetry 
in 4D superspace," JHEP 0203 (2002) 055 |arXiv:hep-th7 0101233|. 

[31] A. Hebecker, "5D super Yang-Mills theory in 4-D superspace, superfield brane oper- 
ators, and applications to orbifold GUTs," Nucl. Phys. B 632 (2002) 101 |arXiv:hep- 
pi/01122 30]. 



33 



[32] N. Seiberg, "Five dimensional SUSY field theories, non-trivial fixed points and string 
dynamics," Phys. Lett. B 388 (1996) 753 |arXiv:liep-th/9608111|; 
K. A. Intriligator, D. R. Morrison and N. Seiberg, "Five-dimensional supersymmetric 
gauge theories and degenerations of Calabi-Yau spaces," Nucl. Phys. B 497 (1997) 56 
tarXiv:hep-th/9702198j . 

A. Hebecker and A. Westphal, "Gauge unification in extra dimensions: Power cor- 



rections vs. higher- dimension operators," Nucl. Phys. B 701 (2004) 273 arXiv:hep- 
[th/0407014j . 



[34] R. Barbieri and A. Strumia, "The 'LEP paradox'," |arXiv:hep-ph /0007265i 

J. A. Casas, J. R. Espinosa and I. Hidalgo, "The MSSM fine tuning problem: A Way 
out," JHEP 0401 (2004) 008 |arXiv:hep-ph/0310137|; 

R. Kitano and Y. Nomura, "A solution to the supersymmetric fine-tuning problem 
within the MSSM," Phys. Lett. B 631 (2005) 58 |ar Xiv:hep-ph/ 0509039|; 
G. F. Giudice and R. Rattazzi, "Living dangerously with low-energy supersymmetry," 
Nucl. Phys. B 757 (2006) 19 (irXiv:hep-ph/0606105] . 

[35] C. Amsler et al. [Particle Data Group], "Review of particle physics," Phys. Lett. B 
667 (2008) 1. 

[36] LEP2 SUSY Working Group [ALEPH, DELPHI, L3 and OPAL experiments], 
| htt p : //lepsusy . web . cern . ch/lepsusy/| 

[37] S. Schael et al. [ALEPH, DELPHI, L3 and OPAL Collaborations, and the LEP Work- 
ing Group on Higgs Boson Searches], "Search for neutral MSSM Higgs bosons at 
LEP," Eur. Phys. J. C 47 (2006) 547 |arX iv:hep -e x/06 02042|. 

[38] G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich and G. Weiglein, "Towards high- 
precision predictions for the MSSM Higgs sector," Eur. Phys. J. C 28 (2003) 133 
|arXiv:hep-ph/0212020j . 

[39] E. Barberio et al. [Heavy Flavor Averaging Group], "Averages of 6— hadron and 
c-hadron Properties at the End of 2007," arXiv:0808.1297| [hep-ex]. 

[40] M. Misiak et al, "Estimate of BR(5 X^i) at C(af)" Phys. Rev. Lett. 98 (2007) 
022002 |arXiv:hep-ph/0609232 |. 

[41] T. Aaltonen et al. [CDF Collaboration], "Search for Bg fJ'^fJ'^ and fJ'^fJ'^ De- 

cays with 2fb-^ oipp Colhsions," Phys. Rev. Lett. 100 (2008) 101802 | arXiv:0712.1708, 
[hep-ex]]. 

[42] J. Dunkley et al. [WMAP Collaboration], "Five- Year Wilkinson Microwave Anisotropy 
Probe (WMAP) Observations: Likelihoods and Parameters from the WMAP data," 
arXiv:0803.0 586 [astro-ph]. 

[43] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, "micrOMEGAs2.0: A program 
to calculate the relic density of dark matter in a generic model," Comput. Phys. 
Commun. 176 (2007) 367 [arXiv:hep-ph/0607059] ; [a^Xi^:0803.2360. [hep-ph]. 



34 



[44] H. Baer, C. h. Chen, F. Paige and X. Tata, "Signals for minimal supergravity at the 
CERN large hadron collider: Multi - jet plus missing energy channel," Phys. Rev. D 
52 (1995) 2746 |arXiv:hep-ph/950327 1|. 

[45] H. Baer, C. h. Chen, F. Paige and X. Tata, "Signals for Minimal Supergravity at 
the C ERN Large H adron Collider II: Multilepton Channels," Phys. Rev. D 53 (1996) 
6241 |arXiv:hep-ph/9512 383|. 

[46] A. Hebecker, J. March-Russell and T. Yanagida, "Higher-dimensional origin of heavy 
sneutrino domination and low-scale leptogenesis," Phys. Lett. B 552 (2003) 229 
(arXiv:hep-ph/0 208249| . 

[47] F. Borzumati and A. Masiero, "Large Muon And Electron Number Violations In 
Supergravity Theories," Phys. Rev. Lett. 57 (1986) 961. 

[48] S. P. Martin and M. T. Vaughn, "Two Loop Renormalization Group Equations For 
Soft Supersymmetry Breaking Couplings," Phys. Rev. D 50 (1994) 2282 [Erratum- 
ibid. D 78 (2008) 039903] |arX iv:hep-ph/9311340 |. 

[49] J. L. Evans, D. E. Morrissey and J. D. Wells, "Higgs boson exempt no-scale supersym- 
metry and its colhder and cosmology implications," Phys. Rev. D 75 (2007) 055017 
|arXiv:hep-ph/0611185] . 

[50] G. A. Blair, W. Porod and P. M. Zerwas, "Reconstructing supersymmetric theories 
at high energy scales," Phys. Rev. D 63 (2001) 017703 |arXiv:hep-ph/0007107 | and 
"The reconstruction of supersymmetric theories at high energy scales." Eur. Phys. J. 



C 27 (2003) 263 arXiv:hep-ph/0210058|. 



35 



